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Abstract 

A combination of analytic and numerical methods has yielded a clear 
understanding of the approach to the singularity in spatially inhomoge- 
neous cosmologies. Strong support is found for the longstanding claim by 
Belinskii, Khalatnikov, and Lifshitz that the collapse is dominated by lo- 
cal Kasner or Mixmaster behavior. The Method of Consistent Potentials 
is used to establish the consistency of asymptotic velocity term dominance 
(AVTD) (local Kasner behavior) in that no terms in Einstein's equations 
will grow exponentially when the VTD solution, obtained by neglecting all 
terms containing spatial derivatives, is substituted into the full equations. 
When the VTD solution is inconsistent, the exponential terms act dy- 
namically as potentials either to drive the system into a consistent AVTD 
regime or to maintain an oscillatory approach to the singularity. 



1 Introduction 

This Chapter, written from the physicist's point of view, will describe an exam- 
ple of synergy between mathematics, physics, and numerical simulation. The 
work reported here was done in collaboration with David Garfinkle (Oakland 
University) , James Isenberg (University of Oregon) , Vincent Moncrief (Yale Uni- 
versity), and Marsha Weaver (Max Planck Institute for Gravitational Physics) 
[ p!9| . Our objective was to explore the approach to the singularity in spatially 
inhomogeneous cosmologies. 

Much earlier (in the 1960's), Belinskii, Khalatnikov, and Lifshitz (BKL) 
[H 1' 0' ^ 1] provided a heuristic analysis of such models suggesting that 
each spatial point would evolve toward the singularity as if it were a spatially 
homogeneous cosmology. We shall call this behavior "locality." In addition to 
the BKL claim of locality, there was also their claim that generically the local 
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approach to the singularity was that of the homogeneous Bianchi Type IX (or 
Mixmaster) model Q . The BKL claims were questioned and it was never 
clear whether or not they were correct ||^. 

To understand these remarks, we begin with the simplest cosmology — the 
spatially homogeneous, isotropic (all points and directions in space are equiva- 
lent) Friedmann-Robertson- Walker models believed to describe our Universe on 
the largest scales Q . Such models are characterized (if one considers evolv- 
ing backward in time) by infinite density, temperature, and tidal gravitational 
force a finite time ago. This is the Big Bang singularity and means that if to 
is the present value of comoving proper time labeling a spacelike hypersurface 
Eq, there exists some labeled by t = such that marks a boundary of the 
spacetime. 

More general, but still spatially homogeneous, cosmologies may be found by 
replacing the isotropic universe expansion rate H = R/R by anisotropic rates 
{Hx, Hy, Hz}- The simplest of these is the Kasner or Bianchi Type I model 
where tJ2^=iHi = 1 = t'^J2i=i^i 1^2 ■ These models are called velocity 



term dominated (VTD) because the only dynamically relevant terms contain 
time derivatives. Spatial derivatives can play a role in spatially homogneous 
cosmologies in the computation of the spatial curvature in a coordinate frame. 
In Bianchi Type I models, the spatial scalar curvature, ^R, is zero. However, it 
is non-zero for the other Bianchi Types. The evolution of these models may be 
understood as dynamics in a configuration space where each point is a cosmology 
at an instant of time. This space will be finite dimensional with directions 
corresponding to each free function needed to describe the space. The space of 
all 3-geometries (3-metrics modulo diffeomorphisms) is called superspace |47|. 
The restriction to homogeneous cosmologies is called minisuperspace (MSS) 1 37 1 . 

The spatial scalar curvature ^R (actually ^R where is the determinant 
of the spatial metric) acts as a potential for the dynamics in MSS. As one 
approaches the singularity, any cosmology in this class behaves as a "free par- 
ticle" — i.e. Kasner model — until a "bounce" off the potential occurs. After 
the bounce, the model behaves as a different Kasner universe. Homogeneous 
cosmologies based on most of the Bianchi Types (in the vacuum case) are asymp- 
totically VTD (AVTD) — there is a final MSS bounce. A fixed set of asymptotic 
Kasner parameters describes the approach to the singularity. We note here that 
topological distinctions among the various Bianchi spaces are dynamically irrel- 
evant. Bianchi Types VIII and IX are not AVTD. The spatial scalar curvature 
forms a closed MSS potential so that there is no last bounce. (Interestingly, 
while this statement is strongly supported by heuristic and numerical evidence, 
there is no rigorous proof. The strongest results so far are those of Rendall 
and Weaver ||44| .) The relationship between one Kasner and the next in subse- 
quent bounces can be expressed as a 4-parameter map which displays sensitivity 
to initial conditions |& 22 , 1^ . Such models are called Mixmaster universes (af- 
ter a food processor) 
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The BKL claim, then, that the approach to the singularity exhibits local 
Mixmaster dynamics (LMD) means the following: if, in some suitable but un- 
specified spacetime foliation, one could write down the general solution to the 
vacuum Einstein equations, the solution at any spatial point as a function of 
time would be arbitrarily close to a Mixmaster universe. The solution looks like 
a Mixmaster universe with spatially dependent Kasner parameters. For this to 
be true, at the very least, the time evolution of the spatial dependence must be 
dynamically negligible compared to the local evolution. The main criticism of 
the BKL claim has been that a specific spacetime slicing is imposed and that 
global issues are not addressed. 

In the work discussed here, a combined analytic and numerical approach pro- 
vides strong support for the BKL claim. The analytic aspect of our methods, 
the Method of Consistent Potentials (MCP) [0, is a "consistency check" for an 
AVTD approach to the singularity. With a suitable choice of variables in some 
spacetime slicing, one may obtain the VTD solution to Einstein's equations by 
solving the system of ODE's which remains if all spatial derivatives are ne- 
glected. In some limit (say r oo), the chosen variables are linear or constant 
in T (presumably with spatially dependent coefficients). The VTD solution is 
then substituted into all the previously neglected terms. If these are exponen- 
tially small as r — + oo, then it is consistent that the model is AVTD. Otherwise, 
one expects that the approach to the singularity is not AVTD. Surprisingly, 
perhaps, the MCP predictions are borne out by detailed numerical simulations 



of the full Einstein equations for these models ]19|, |16|, |l4 |17|, [46[. 

In recent years, it has become possible to prove that the behavior of certain 
simple spatially inhomogeneous cosmologies is AVTD. An example is the po- 
larized Gowdy model consisting of a single polarization of gravitational waves 
with symmetry in the plane perpendicular to the direction of propagation. The 
wave amplitude satisfies a linear equation while the rest of the spacetime may 
be constructed if the wave amplitude is known. The addition of the second 
polarization completely changes the character of the solution. The two wave 
amplitudes satisfy nonlinearly coupled wave equations. Numerical solutions, as 
interpreted through the MCP, for thses models have allowed sufficient under- 
standing of the solutions so that rigorous study is possible. It can now be proved 
that solutions with the numerically observed behavior exist [Q. This is not a 
redundant effort because the proof overcomes objections concerning numerical 
artifacts or limitations on the ability to explore initial data space. Furthermore, 
the numerical studies cannot be carried to r = oo. 

In Section II, we shall present an introduction to homogeneous cosmologies 
including a detailed example of the use of the MCP. In Section II, we shall 
discuss Gowdy models — showing how some are AVTD, interesting features of 
the approach to the singularity, and how adding a magnetic field yields LMD 
[ p^ [l^ , In Section IV, we discuss progress to date on the most complicated 
cosmologies we have studied — those with a single spatial C/(l) symmetry (which 
I propose to call Moncricf models) |l^ — and show the relation to the 
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homogeneous Mixmaster model. Discussion will be given in Section V. 



2 Spatially homogeneous cosmologies 

Bianchi's classification of homogeneous 3-spaces implies that they can be 
described by metric 1-forms, fT% such that (for i,j — 1,2,3) 

da' = q^CT^ A a'' (1) 

where C*^. are the structure constants of a group (e.g. SU{2) for Bianchi Type 
IX). A homogeneous spacetime then has the metric 

ds^ = -N^dr^ + e^'' (e^^) a' a' (2) 

where 

A, = diag (-2/3+, p+ + V3I3- , 13+ - . (3) 

The lapse N{t) allows a change of time coordinate. The spacetime foliation 
is chosen to make the homogeneity manifest. We shall follow BKL and take 
jY _ g3n_ rpjj^g chosen variables, 51 (r) and /3±(t) respectively, measure the 
volume and anisotropy (two orthogonal shears) of the space at a given time. 
MSS is the 3-dimensional space with axes (fi, /3_). Each point in MSS is a 
homogeneous space. A trajectory in MSS is a cosmological spacetime. 

To understand the dynamics in MSS, it is most convenient to use the Hamil- 
tonian formulation of Einstein's equations pioneered by Arnowitt, Deser, and 
Misner (ADM) 0, H, When speciahzed to homo geneous spaces, Einstein's 
equations for Class A Bianchi types (the only ones we shall consider here) may 
be obtained by variation of the "superhamiltonian," NH, where H = is the 
Hamiltonian constraint. For the Kasner and Bianchi Type IX Mixmaster models 
we find 

2H = -pI +p\+ pi + y (/3± , n) (4) 

where 

, , f , Bianchi Type I 

[Vix , Bianchi Type IX 

for 

-2e4^+ - 2e-2(^++^'3-) _ 2e-2(/3+-v^/3-)] 

while {pn,P±) are the momenta canonically conjugate to {^,P±)- The VTD 
solution is just the Kasner solution 

n = Qq- pnT , 
P± = Pl+P±r (7) 
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with (from H = Q) 

pI=p\+pI. (8) 

To apply the MCP to the Bianchi Type IX model, substitute the Kasner solution 
into (|) to give (as r oo) 

T/lim _ 4|pr2|T(-l+2cose) I _^4|pQ|r(-l-cos0+\/3sine) 

Vjx — e + e 

_|_g4|pn|T(-l-cose-y3sin6l) ^g'j 

where cos 9 — —p^/pn and sin 9 — —p^/pn. The terms indicated by . . . are 
almost always negligible j2| |ll] and we shall ignore them here. Clearly, except 
for the set of measure zero 9 = (0, ±27r/3), one of the terms on the right hand 
side of (H) will grow. This leads to the MCP prediction that generic Bianchi 
Type IX models are not AVTD. This, of course, is consistent with the belief 
that there is no last bounce in these models. 

Asymptotically, the Bianchi Type IX model solution may be described as a 
presumably infinite sequence of Kasner parameters. It was shown long ago by 
BKL Q that the Kasner exponents {pi}, such that the anisotropic scale factor 
Hi = t^'* , may be expressed in terms of a single parameter 1 < u < oo via 

Pi = —u/ [u^ + u + 1) , 

P2 = (U+ l)/(u2 +M+ 1) , 

Pi = u{u+l)/{v? + U+1) (10) 

where 

for Un the u-parameter of the rtth Kasner epoch. This map contains the infor- 
mation in the Gauss map 

UN+l = \ T (12) 

UN - [un\ 

where [ ] denotes integer part. The Gauss map is known to be chaotic Q]. To 
arbitrary accuracy, as long as the terms indicated by ... in (^) are negligible, an 
approximate evolution may be constructed using the 4-parameter map ^ ^ 
as a sequence of Kasner models with the rules to go from one Kasner model to 
the next. From this analysis, one can approximate the duration of each Kasner 
epoch (in r or \^\). This quantity depends on the values of the parameters but 
increases roughly exponentially in r toward the singularity. 

If the BKL claim is correct, asymptotically each spatial point would be 
characterized by a spatially dependent u{x)^ obeying ( pi] ) to ever better accuracy 
in the approach to the singularity. But the Kasner epoch duration dependence 
on u{x) means that bounces occur at different times at different spatial points, 
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leading eventually to a complicated and arguably fractal spatial dependence of 
the metric [g |§. 

For later reference, we note here that the addition of magnetic fields to 
AVTD Bianchi Type vacuum cosmologies can "close" the MSS potential to yield 
a Mixmaster model [^. One example is the magnetic Bianchi VIq cosmology 
with H, |l| 

^g40+4/3+-4v^/3_ _^ ^ ^ ^ (13) 

where . . . indicates terms different from those in (^) but still dynamically irrel- 
evant almost always. 



3 Gowdy cosmologies and their generalizations 

The simplest solutions representing spatially inhomogeneous cosmological space- 
times are the polarized Gowdy models onT^ xR ^ described by the metric 

where A and P are functions of < 6' < 27r and time, r, only. Here P satisfies 
a linear wave equation 

P,rr p,eg = (15) 

with explicit solution 

oo 

P = Po + aoT + ^ Q;„Zo(ne~^) cos(n6' -I- 0„) (16) 

n=l 

where Z^{x) is a zero-order Bessel function. Given P, A may be constructed 
from the Hamiltonian and momentum constraints written respectively as 

A,.= -(P,.)'-e-2-(P,,)2, (17) 

\,e - -2P,r P,e ■ (18) 

This model may be obtained as an inhomogeneous generalization of a Kasner 
solution. The approach to the singularity was shown heuristically early on [ p^ 
and later rigorously to be AVTD. It was also proven that the solution may 
be extended to the strong curvature singularity at r = oo and only becomes 
singular in that limit |2^ . 

The Gowdy models become more interesting when the second gravitational 
wave polarization is added [l^ . The metric becomes 

ds" = e^/2(-e-3^/2^r2-He^/2^02) 

-fe-^[e^dCT2 ^ 2QePdadS + [Q^ + e"^) dS"^] (19) 
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where P and Q satisfy the nonhnearly coupled equations 

P,rr -e-^^P,ee -e^''{Q,l -e'^^Q,l ) - 0, (20) 

Q,rr -e-^^Q,ee +2P,r Q,r ~2e-^^P,e Q,e = (21) 
while the constraints become 

A,. = -P,l -e'^Q,l -e-'^P,l +e'^Qi ), (22) 

X,e=-2{P,rP,e+e^^Q,rQ,e). (23) 

This model is still an inhomogeneous generalization of the Kasner solution but 
with a rotation of the principle axes of the 3-space between one spatial point 
and another. 

To apply the MCP, it is convenient to use the Hamiltonian density 

2H=nl+ e-^^nl + er^^ (Pj + e^^Q J ) (24) 

whose variation yields the wave equations (^^ and (pTf). The first step is to 
solve (20)-(|3|) neglecting all terms containing spatial derivatives. In the limit 
as r cx), we find the VTD solution ||l6|, ^ 

P = VT , Q = Qo , A = -wV (25) 



where v depends on 6 but not on t. Now substitute ( p5| ) into (24). We see that 

V, = \nle-'^ (26) 
can be exponentially small if d > while 

V2 = lQ3e'^~'^ (27) 

can be exponentially small if w < 1. Thus the MCP predicts that the approach 
to the singularity in these models will be AVTD (almost everywhere [^) if 
0<v <1. 

The question then arises as to what happens if w < or w > 1 in the 
initial data. Prior to the asymptotic regime, spatial derivatives need not be 
dynamically irrelevant and a variety of nonlinear effects may occur . At some 
value of T at a given = Oi, the true Hamiltonian is excellently approximated 
by either 

Hi = + \^le-^^ (28) 

or 

H.^\-^1 + \q.I^''~'^. (29) 
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In the first case, there will be a bounce off Vi causing np (and thus v{Oi)) to 
change sign. But Vi will cause a bounce only when it is exponentially growing 
— i.e. when v{9i) is initially < 0. The bounce will then yield the new value 
v'{9i) > 0. In the second case, v{9i) — 1 changes sign at the bounce off V2. 
Here the bounce occurs only for v{9i) > 1 so that V2 grows exponentially. If 
v'{9,) - 1 = -iv{9,) - 1), then v'{9,) =2 - v{9,). If v'{9,) < 0, there will 
be another bounce off Vi and so on. No additional bounces will occur once 
< v{9i) < 1. Thus the MCP predicts a sequence of bounces off Vi and V2 at 
spatial points where v{9) < or v{9) > 1 until < v{9) < I at which point the 
true solution approaches a VTD solution at these points. 

There is a set of measure zero at which the coefficient of Vi (i.e. ttq) or 
V2 (i.e. Q,e) vanishes. At these spatial points, the above mechanism fails to 
operate so that v{9) > 1 or v{9) < is allowed Q. If, at 6*0, Q,e = 0, then, for 
9 « 9o, Q,e ~ so that v{9) > 1 can persist for a long time. This creates "spiky" 
features in P in the vicinity of one of these non-generic points. Similarly, P,^ < 
can persist if ttq « 0. This in turn leads to an apparent "discontinuity" in Q 
since Q,t = ttq e~^^ grows exponentially in the positive or negative direction 
about the region where ttq goes through zero with P,,- < 0. 

All the properties described here in terms of the MCP have been observed 
in numerical simulations of the full Einstein equations (|2^)-(|3|) This 
provides strong support for the BKL claim that each spatial point evolves toward 
the singularity as a separate universe. So far, these models are too restrictive 
to admit LMD. Once the behavior of these models became apparent in the 
numerical simulations, Kichenassamy and Rcndall | |33[ | were able to apply the 
former's methods to prove that solutions with this property exist. Although the 
numerical simulations provide information about the nature of generic solutions, 
they cannot explore all initial data space nor can they evolve to t = 00. The 
MCP provides a framework to interpret observed waveforms and thus to identify 
generic behavior independent of initial conditions. Similarly, it allows one to 
argue that nothing qualitatively new will happen if the simulations are followed 
longer. While the studies in |Q do not prove the solutions are AVTD (in 
contrast to the result js^ for polarized Gowdy models), they do provide some 
rigorous verification of the numerical results. 

To overcome the limitations of the vacuum Gowdy models, consider the 
generalization of the metric (|l^) to represent the inhomogeneous version of 
the magnetic Bianchi VIo model |46[ It is necessary to change the spatial 
topology from to the solv-twisted torus with boundary conditions ^ 
X{9 + 2Tr,T) = X{9,t), P{9 + 2tt,t) = P{9,T) + 27ra, Q{9 + 2ti,t) = Q(9 ,t) e''^'^"- 
where a is a constant. To incorporate the coupling to the magnetic field, the 
metric must be generalized and becomes 

ds" = e^/2(-e-3-/2dr2+e(A<+-)/2d^?2) 

+e-^ {e^da'^ + 2Qe^dadS + (Q^ + e-^) d6^] (30) 
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which yields the Hamihonian density 

TT 1/2,2 -2P I -2rr)2 , 2(P-r)^2A 

+47rAe(^+^'/2 ^2 (31) 

where _B is a constant related to the magnitude of the magnetic field. The 
constraint Tr\ = ^e'^/^ has been used to obtain this expression The 
VTD solution is the same as (^5|) for the vacuum case if we make the substitution 

vie)^P,r^^^ (32) 

where ttx also depends on 9 but not on r. The new definitions of Vi and V2 
require only division by 2'!T\. All the arguments regarding Vi and V2 go through 
as before. However, if < v{9) < 1, the new potential in (^l|), 

V3 = 47rAe(^+")/' , (33) 

will grow exponentially. Thus the dynamics will be dominated by one of Hi 
(from (H)), H2 (from (|2£ 



i/3 ^ + 47rAe(^+")/' . (34) 

Since no value of v{0) is consistent with Vi, V2, and V3 all exponentially small, 
the MCP predicts that the magnetic Gowdy model is not AVTD. Since this 
magnetic Gowdy model was obtained from magnetic Bianchi Type VIo, the 
observed departure from AVTD behavior should indicate LMD. Numerical sim- 
ulations support the MCP predictions to the extent that they can be followed 
So far, nothing rigorous is known about the solutions. 

4 ^(1) symmetric cosmologies (Moncrief mod- 
els) 

Some years ago, Moncrief used an ADM-like approach to write Einstein's equa- 
tions for a cosmology with a single spatial U{1) symmetry |Q. He considered 
Einstein's equations projected along and normal to the Killing direction. It is 
convenient to use a metric of the form 

ds^ = j-iVdr^ + gabidx" + N''dT){dx^ + TV^dr) j 

+e^'^ [dx^ + Padx" + Podr) ^ (35) 

where x^ is the Killing direction, x^ — u and x'^ — v are the other spatial 
coordinates, while r is the time. The variables are ip^ where e"^ is the norm of the 
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Killing field and /3" for a = 1, 2 are the "twists." The conformal lapse N (chosen 
to be iV = e^), shifts N"" (chosen to vanish), and /3o become Lagrange multipliers 
of the constraints. The spatial U{1) symmetry allows the last constraint, e°,a = 
0, to be satisfied identically through the introduction of the twist potential lu 
such that e° — e"'''uj,b where e"^ — —e^°-. Thel 2-metric in the u-v plane is 
parametrized as 

9ab = e^eab (36) 

with 

_ 1 /e2^ + 6-2^(1 + a;)2 e^' + e-^'{x^ - 1)\ . . 

eafc- 2 l^e2. ^g-22(^2_i) g2. ^g-2.(i_^)2 j ■ l^Tj 

If the variables {if, lu, A, x, z} are joined by their canonically conjugate momenta 
{p,r,pA,px,Pz}, it is found that Einstein's equations may be obtained by vari- 
ation of the Hamiltonian density ||l^, |l^ 

Tj -i2|-i42 2|-i2il4M2 ^2 

+ (e V") ,ab - (e^e"") ,a A,b [(e-^-) ,„ x,, - (g-^-) x,u ] 
+2e^e'^V,a +^e^e-^^e'''u;,a uj,b (38) 

where iJ = is the Hamiltonian constraint. Since the structure of this Hamil- 
tonian is reminiscent of that for the Gowdy models (HJ) and (|3ll), one may 
consider application of the MCP. The VTD solution is 

Z = -VzT, X = Xo, p^ = -4Vz, P^=pI, ip=-V^T, 

Lo = Wo, p = —iv^, r^r°, A = Aq - wat, pa = wa (39) 

where Vz, v^, Xq, Wq, , Ag, and v\ > are functions of u and v but 
independent of r. (The sign of ua is fixed to ensure collapse.) Although quite 
complicated, the nonlinear terms in ( ^ ) fall into four classes: 

(1) First, the term 

Vz = \ple''^ (40) 

grows exponentially only if < 0; 

(2) Second, we find that 

^1 = ^r^e^'^ (41) 

grows exponentially only if < 0. 

If we then assume Vz > > 0, the inverse conformal 2-metric is 

dominated by e"^^ which grows exponentially. 

(3) It is then seen that all but one term in ( |38| ) containing spatial derivatives 
has the factor 

F = e^~'^' . (42) 
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(4) The remaining exponential term is 



The VTD form of the Hamiltonian constraint is 



i?VTD ^-\pI + + \p'r- = - ^^'A + K + = • (44) 



For > to ensure collapse, > 0, and > 0, (^44| ) means that 
e^~2^ will always be exponentially small since i;a > 21;^ while V2 will always 
grow exponentially since da < 2^2 + ^v^. This means that > 0, > 
is inconsistent and a bounce off V2 will always occur. The "bounce rules" — 
equivalent to conservation of momentum — are pO| 



= -v^ + -UA - , 



1 

—1 
2 

1 1 



-v^ + -VA + -jV^ 

1 

—1 
2 



v'a -2wip + -UA - t^z (45) 



indicating that changes sign at the bounce. Occasionally, if Vz > ^v^p — 
i^A, Vz can also change sign. If either Vz or v^p is < 0, a bounce off Vz or 
Vi will subsequently occur to reverse the sign. Thus the MCP predicts an 
infinite sequence of bounces (almost) everywhere. (The behavior near non- 
generic spatial points has not yet been studied for these models.) 

The condition r = w = is preserved by Einstein's equations and yields po- 
larized Moncrief models. In these models, both Vi and V2 are identically zero. 
If < initially at some spatial point, the MCP predicts that a bounce off Vz 
will occur to yield z increasingly large and negative. Since the surviving expo- 
nential terms with factor F are exponentially small under these circumstances, 
the MCP predicts that polarized Moncrief models should be AVTD. Numerical 
simulations for these models do not have the problematic numerical aspects of 
generic Moncrief models and yield strong evidence to support the MCP predic- 
tion 0. The methods of Kichenassamy can be used to prove that these 
models are AVTD 

The constraints in the Moncrief models are non-trivial and must be solved 
initially. In published work, we have used the following algebraic prescription 
to solve the constraints jlj]: Solve the momentum constraints (see [^0[ |l^) by 
imposing Px = Pz — — — and pA = ce^. For c sufficiently large, the 
Hamiltonian constraint may be solved for p or r. Thus x, z. A, and r or p may 
be freely specified. 

Numerical simulations of these models are less straightforward than in the 
Gowdy case. In order to be able to follow the simulations toward the singularity, 
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it appears to be necessary to prevent the development of Gowdy-like spiky 
features by introducing data smoothing — every variable is replaced by 
a 6th order accurate average of itself over a grid of nearby spatial points. In 
addition, the crucial and explicit role of -ffvTD in the dynamics means that it 
is essential to preserve the Hamiltonian constraint. This is done by solving the 
Hamiltonian constraint, H (from (^8|)) — 0, forpA at every time step. We believe 
that both these numerical procedures not only preserve a stable evolution and 
but also preserve the qualitative features — especially the evolution in time at 
a fixed point in space — of the actual solution 

Once again the numerical simulations support the MCP predictions. This is 
perhaps the most remarkable case since, for the MCP to yield valid predictions, 
the time scale for the change in the coefficients of the exponentials in Vi , K: , 
and V2 must be significantly longer than that for the exponential itself. 

Aside from the need for improved numerics, the remaining major open ques- 
tion is whether the oscillations in Lp observed in Moncrief models correspond 
to LMD. This issue may be addressed by writing a homogeneous Mixmaster 
model (on S'^) in the Moncrief (on S^) variables — i.e. identify the relation 
between the variables in the metric ( ]35| ) and those in the metric (|^) by writing 
the metric 1-forms cr* in a coordinate basis as = cos (f>d6 + sin 6 ain 4> dip ^ 

= sin (j)d6 — sin cos (j) dip, — dcp + cos Odip |3^. It is then straightforward 
to identify the cyclic coordinate ip with the U{1) symmetry direction x^. Details 
of the transformation between Moncrief and MSS variables are given elsewhere 
[|8[|0). 

For our purposes here, it is important to consider the Mixmaster models in 
terms of the logarithmic scale factors 

a = n-2(3+, 

C = r! + /3+ + V3/3_ , 

7 = n + l3+^V3(3- (46) 

— i.e. e" is the 1-direction scale factor, etc. In a regime such that a > C > 7, 
then 



ip = a + —\n 



sin^ e sin^ 6 + e^^-^-") cos^ 6 sin^ 9 + e^^^"") cos^ ( 



(47) 



This expression is exact but is written to express the dominance of the largest 
logarithmic scale factor (LSF). The result for any other ordering of the scale 
factors is obtained by the appropriate permutation. The identification of with 
the largest LSF creates an oscillatory behavior in ip which may be understood 
from the behavior of the LSF's in Mixmaster dynamics. For a > ^ > 7, a 
increases while C decreases until there is a bounce. At this bounce, a and ^ 
(where ' = d/dr) change signs so that a decreases and C, increases. The Moncrief 
variable 93 follows a as long as a > C — through a sign change of ip^^ at the 
bounce. When a — there is another sign change in Lp^^ as it starts to track 
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C,. While this is going on, 7 decreases monotonically. At the bounce, 7 becomes 
less negative. Eventually (when 1 < u < 2 in (|Tl|)), 7 will change sign. 

The Moncrief variable z is written in terms of logarithmic scale factors as 



sin2 6'|sin<?!)|(l + ...) 



1 - cos 29 



(48) 



where terms indicated by . . . are exponentially small if a > C > 7- We see 
that z is large and negative unless 7 ~ C — i-^- 7 > so that 7 is increasing. 
This marks the end of the BKL era. Observation in numerical simulations of 
an increase in z would provide a qualitative signature for LMD which would be 
independent of numerical uncertainties. 

The usual Mixmaster bounces show up in the Moncrief models as bounces 
of (yC off Vi. The "extra" bounces in Lp are the bounces off V2 which is not 
exponentially small only when a « C (for our original ordering of the LSFs). 
The bounce in z to halt its descent to large negative values also occurs off V2 
and can be understood from the V2 bounce laws (^5|). The increase in z will 
eventually be halted by a bounce off ■ 



5 Conclusions 

Despite the complexity of Einstein's equations, strong statements may be made 
about the approach to the singularity. To make these statements, we have 
obtained a synergy among heuristic, numerical, and mathematical analyses. 
The MCP has been applied to spatially inhomogeneous cosmologies with as few 
as one spatial Killing vector field to predict whether or not the approach to the 
singularity is AVTD. The observed numerical behavior has already been used in 
the vacuum Gowdy and polarized Moncrief cases to guide the proof of theorems. 

The original BKL picture remained difficult to understand because BKL 
used more of the VTD solution than was really necessary. Our consideration of 
the consistency of linear asymptotic behavior in the MCP clarified the picture 
considerably. It is now possible to determine quickly, given the appropriate 
choice of variables, how a collapsing model should behave. 

This brings us to the question of whether our results are dependent on the 
choice of variables and the spacetime slicing. As far as the former goes, we 
now have two completely different descriptions of the Mixmaster universe — 
one in terms of the MSS variables (3± , and the other in terms of the Moncrief 
variables <p, z, etc. In both cases, the VTD solutions for the variables are linear 
or constant in r. (The same spacetime slicing is used in both descriptions.) 
In both cases, there are exponential potentials which cannot be kept small. 
However, as we have seen, the meaning of the variables and potentials is not 
the same in both cases. When applied to the inhomogeneous models, we may 
conjecture that, even within a given spacetime foliation, there may be many 
choices of variables which can yield an MCP prediction and that the prediction 
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— AVTD or not — will be independent of the choice. Certainly, one should be 
able to construct analogs of the MSS variables for the Moncrief models. 

One might further conjecture that all reasonable foliations will converge 
asymptotically in the sense that, e.g. in a solution which is AVTD, in some 
foliation, it will be possible to choose a time coordinate such that there are 
local Kasner variables linear in that time coordinate. Some analytic work on 
polarized Moncrief models and Gowdy models supports this conjecture. 

We emphasize here the synergy between the MCP, the numerical simulations, 
and the mathematics. To have only the simulations would suffer the limitations 
of restricted initial data and limited duration. The MCP provides an organizing 
principle to allow conjectures to be made about generic behavior based on the 
validity of the MCP as applied to the simulations. Conjectures about the generic 
behavior then form the basis for mathematical analysis. 

So far, we have explored inhomogeneous cosmologies with one and two 
Killing fields. The observed local behavior in all cases suggests that the bound- 
ary conditions are not important and that the results are applicable in any 
collapse — e.g. in asymptotically flat spacetimes. Some support for this exists 
in the observations of typical Gowdy behavior in S'^ x rather than spatial 
topology [ p5[ and in the same generic Mixmaster behavior in magnetic Gowdy 
models on the solv- twisted torus or |l3|. 

One remaining open question in the Moncrief models is the extent to which 
the failure to be AVTD indicates LMD. Writing a Bianchi Type IX Mixmaster 
cosmology as a Moncrief model provides a qualitative signature for LMD in the 
rise of z. This may have been observed |Q. So far, numerical difficulties make 
it impossible to conclude that every feature in the asymptotic regime is LMD. 
Further studies are in progress. 

What about the zero Killing field case? The local behavior may allow this 
problem to be tractable numerically since high spatial resolution would not be 
required. There is no obstruction to use of the MCP in this case. In fact, 
Moncrief has obtained the Hamiltonian for this sytem starting from the metric 
( ^5|) but assuming dependence on as well as u, v, and r. The transformation 
to the twist potential (r, u) now cannot be made yielding the expected six 
degrees of freedom {ip, x, z, A}. 

Every term present in the Moncrief model is present in the more general 
case so that LMD should be possible. It is not yet clear to what extent there 
are additional terms which may signal a richer asymptotic behavior or whether 
the BKL picture will still hold in the most general case. 
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